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Abstract 

l/^ [ A Bargmann symmetry constraint is proposed for the Lax pairs and the adjoint Lax 

, pairs of the Dirac systems. It is shown that the spatial part of the nonhnearized Lax pairs 

and adjoint Lax pairs is a finite dimensional Liouville integrable Hamiltonian system and 
■ that under the control of the spatial part, the time parts of the nonhnearized Lax pairs 

Q ' and adjoint Lax pairs are interpreted as a hierarchy of commutative, finite dimensional 

' Liouville integrable Hamiltonian systems whose Hamiltonian functions consist of a series of 

integrals of motion for the spatial part. Moreover an involutive representation of solutions 
fS| ■ of the Dirac systems exhibits their integrability by quadratures. This kind of symmetry 



o 
o 



X 



constraint procedure involving the spectral problem and the adjoint spectral problem is 



' referred to as a binary nonlinearization technique like a binary Darboux transformation. 

' Key Words: Zero Curvature Representation, Nonlinerization Method, Liouville In- 

tegrable system, Soliton Hierarchy 

1. Introduction 

O ■ Symmetry constraints become prominent in recent few years due to the important roles they 

play in soliton theory. For 1 + 1 dimensional soliton hierarchies, a very successful symmetry con- 
straint method is the nonlinearization technique proposed by Cao^^'^l and Cao and Gengl'^'^l . A large 
class of finite dimensional Liouville integrable Hamiltonian systems is thus generated which are con- 
nected with soliton hierarchies^^"^] and the nonlinearization technique is systematically extended 
by Zeng and Li^^'^l for quite a few soliton hierarchies. At the same time, a symmetry constraint 
procedure for bi-Hamiltonian soliton hierarchies is proposed by Antonowicz and Wojciechowskil^l 
and bi-Hamiltonian structures for the constrained systems can be worked out^^'^^l. By observing 
that stationary soliton fiows may be interpreted as finite dimensional Hamiltonian systems^^^l on 
introducing the so-called Jacobi-Ostrogradsky coordinates!^^! , a natural generalization of nonlin- 
earization technique to higher order symmetry constraints is made by Zengl^'^1 for the KdV and 
Boussinesq hierarchies. 

For 1 + 2 dimensional soliton hierarchies, a sort of symmetry constraints similar to ones in 
the nonlinearization technique has also been presented by several authors^^^], although relatively 
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little is known about the general construction at present. Some interrelations of 1 + 2 dimensional 
soliton hierarchies to 1 + 1 dimensional soliton ones and various structures of constrained systems 
have vigorously been exposed^^^l . The theory progress along the above two lines has been making 
and has aroused increasing interest in recent times. 

This paper is devoted to the symmetry constraints in the first line. Based upon the idea of 
nonlinear ization technique, we would like to propose a binary nonlinear ization technique for soliton 
hierarchies possessing Lax pairs. Wc shall express, in terms of the eigenfunctions and the adjoint 
eigenfunctions, the variational derivative of the spectral parameter with respect to the potential, 
and thus the binary nonlinearization technique involves two sets of dependent variables, which is 
different from the nonlinearization technique. By considering the constraint problem of the Lax 
pairs and the adjoint Lax pairs, we shall exhibit the idea of binary nonlinearization technique and 
provide new examples for finite dimensional Liouville integrable Hamiltonian systems. In Section 
2, we shall give the concrete construction of the Dirae systems and their Hamiltonian structures. 
Furthermore we want to establish some properties which will be used later. In Section 3, we shall 
exhibit systematically the binary nonlinearization procedure for the Lax pairs and the adjoint Lax 
pairs of the Dirac systems. 

2. The Dirac systems and their Hamiltonian structures 
We consider the Dirac spectral probleml^^^ 

cl>, = Ucf> = Uiu,X)cP, <^=(^^)' ^=(r)' (^-^^^ 
where the spectral operator is as follows 

U = U{u, A) = rai + Xia2 + qas = ^ ^ , (2.1b) 

and the aj, 1 < j < 3, are 2x2 Pauli matrices. To derive Dirac systems associated with (2.1), we 
first solve the adjoint representation equation (see Ref. [17]) Vx = [U, V] of cpx = U(p. Set 

V = aai+ bia2 + c^s = ^ ^ ^ "^^^ j . (2.2) 

Noting that 

[U, V] = (-2Ac + 2qb)ai + (-2rc + 2qa)ia2 + (2Aa - 2rb)a3, 

we see that the adjoint representation equation T4 = [U, V] becomes 

' = — 2Ac + 2qb, 
bx = —2rc + 2qa, 
^ Cx = 2Xa — 2rb, 
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(2.3) 



which is equivalent to 

' ao = Co = 0, 6ox = 0, 
ttix = -2cj+i + 2qbi, 
hx = -'^rci + 2qai, i > 0, 
^. Ci^ = 2ai+i - 2rbi, 

on setting a = Yli>o Oi^~^-, b = J2i>o bi^~^-, c = Ylii>o Ci A~*. We choose ho = —1, and assume that 
o-i\u=o = bi\u=o = Ci|u=o = 0) ^ > 1 (or equivalently select constants of integration to be zero). At 
this point, the recursion relation (2.3) uniquely determines a series of differential functions with 
respect to u. For instance, we have 



ai = 


-r, ci 




1 


a2 = 






1 


as = 


■^r^x - 



= -q, hi = 0; 

C2 = -^r^:, 62 = 
1 2„ 1 3 , 



-^2 ^2 

■2^ -2^' 

1 _ 1 3 

- ^Qxx 2^ 



2" 2 ' 
From (1/2)^ = we see that (^trF^ 

have + — 6^ = 1. Further we obtain 

^ n — l 

bn ~ ^ '] {bibn—i Oiflra— i CjC^— i)) Tl > 2. 



(c2 + a2 - 62), = 0. Thus by atiV^)\u=o = h we 



By the mathematical induction, it follows from (2.3) and the above equality that ai,bi,Ci are all 
differential polynomial functions of u. The compatibility conditions of Lax pairs 



(t>^ = U(f>, = y(-)<^, = {X^V)+, n > 0, 



(2.4) 



where the symbol + stands for the selection of the polynomial part of A, engender a hierarchy of 
the Dirac systems 

-2ari_|_i 



2a 



■n+l 



n > 0. 



(2.5) 



The first nonlinear Dirac system in the hierarchy (2.5) reads as 



9*2 



-^r^^ + q'^r + r^, 

-Axx -q^ - qr^- 



This system is different from the coupled nonlinear Schrodinger systems in AKNS hierarchy because 
it contains the cubic terms q^,r^. 

In what follows, we want to give the Hamiltonian structures of the Dirac systems (2.5) by 
means of the trace identity proposed in Ref. [18]. To this end, we need the following quantities 
which are easy to prove: 



OA oq or 



(2.6) 



where < A,B >= tr(AB). Now applying the trace identityl^^l 

S ^.dU , „ a „ dU „ dU 



we obtain at once 



Equating the coefficients of A ^ on two sides of the above equality, we have 

-^{-bn+i) = {7 - n)icn,an, f , n > 0. 

By taking simply n = 1, we find that the constant 7 = 0. In this way we obtain an important 
equality 

A/, 

-- n{cn,an,) , n > 0. (2.7) 



6u 

In addition, for ra > 1 we have 

2cn+i ) \l Q )\2an+i) \l Q )\\d-2rd-^r 2rd-^q J \2an 
Therefore the hierarchy (2.5) may be cast into the following Hamiltonian form 



where the Hamiltonian operator J, the recursive operator L and the Hamiltonian functions are 
determined by 

By a direct calculation similar to Ref. [19], we can get 

y/™) - y/^) + = (F('"))'[i^„] - (yW)'[i^„] + yW] = o, 

which implies the commutability of the flows of (2.8). Therefore each system in the hierarchy (2.8) 
has an infinite number of symmetries {Km}m=o- Besides, we can directly verify 

Vt„ = [V^''KV], n>0, (2.10) 

when ut„ = K^, i.e. Ut„ - V^"^ + [U, = 0, n > 0. In fact, we easily find that Vt„ - F] 
satisfies the adjoint representation equation of (f>x = U(f) and that T4„ — F] vanishes at u = 0. 

Thus (2.10) holds forn > because the adjoint representation equation T4 = [U, V] has uniqueness, 
namely if 14 = [U, V] and V vanishes at u = 0, then V itself vanishes. 
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3. Binary nonlinearization 



Let us begin with the binary nonUnearization of the Lax pairs and adjoint Lax pairs of the 
Dirac systems. Associated with Lax pairs (2.4), the Dirac systems have the adjoint Lax pairs 

= -U^i; = -U^{u, (3.1a) 
{ = -(yW)^V = -{V^^^)'^{u, A)V, (3.16) 
where T means the transpose of matrix and = (V'l) '02)^- It follows from (p^ = Ucf), ipx = —U'^i^-, 
that 

6\ 1 

-T- = —{(pllpl- (I)2tp2,(pl'ip2 + 4'2'ipl)^, E = {(pltp2 - (p2'ipl)dx. (3.2) 

ou E 

When the zero boundary conditions: lim|^|^_|_oo (f) = lim^x^_,_^_^ ip = 0, hold, we have 



L— = A—, or L J— = AJ— , 

ou ou Ou Ou 



(3.3) 



where L is defined as in (2.9), and L* is the adjoint operator of L, which is a hereditary symmetry. 
For a general spectral problem, this kind of property (3.3) has been discussed in Ref. [20]. 

Now introducing distinct eigenvalues Ai, A2, • • • , Xn, we obtain the following two spatial and 
time systems 



4>l3 



U{u, Xj 



^•^1 i = 1 2 ■ ■ ■ N 

&2j. 



-C/^(u,A,) 



i = 1,2, •••,7V; 



^■^3 / t 



F(-)(^,A,)(^^jV i = l,2,...,iV, 



= -(F("))^(u,A,)( 



, j = l,2,---,iV. 



(3.4a) 

(3.46) 
(3.5a) 

(3.56) 



Because Ut^ - + [U, F^*^)] = if and only if {-U^)t„ - (-(V'^"))^)^ + [-U^^\ -(F^"))^] = 0, 
the compatibility condition of (3.4) and (3.5) is still the nth Dirac systems Ut„ = Kn- Let us take 
the Bargmann constraint (requiring the G-vector field to be a linear function but not a differential 
function with respect to u) for the Dirac systems 



JGo = jY^^2Ej-^, Ej = J {(t>iji'2j - hji'ij) dx. 
of the Dirac systems. The constraint (3.6) allows us to impose that 



(3.6) 



This kind of constraints is, in fact, symmetry constraints because the J-r^ are common symmetries 

ou 



^^^^ = (< ^1,^1 > - < $2,^2 >,< ^1,^2 > + < ^2,*l >)^, 



which impUes that 

j q = - <^i,^i> + <^2,^2>, (3.7a) 
I r = - < $1, *2 > - < ^2, *i > • (3.76) 

Here $j = {4>ii, ■ ■ ■ , (l)iNV , = (V'ii, • • • , V'iiv)'^, i = 1,2, < y,z >= YJj^iVjZj, y,z G . 
Note that (3.7) has two particular properties: including both the eigenfunctions ^1,^2 and the 
adjoint eigenfunctions ^'i, ^'2, and nonlinearity with respect to Therefore we refer to (3.7) 

as a binary nonlinear constraint. In general, nonlinear constraints only involve eigenfunctions 
of spectral problems associated with integrable systems (see Refs. [1-8]). We denote by A the 
expression of A under the constraint (3.7). The property (3.3) ensures that 

/ < *i > - < A"-i$2, *2 > ' 

\ < *2 > + < ^""^$2, *1 > , 

and that from (2.3), 

K = d-\2qan - 2f c„) = < ^""^$1, *2 > - < ^""^$2, ^'i >, n > 1. (3.8b) 

Here bi 7^ 0. But = [U, V] still holds. By substituting (3.7) into the Lax pairs and the adjoint 
Lax pairs: (3.4) and (3.5), we acquire the nonlinearized Lax pairs and adjoint Lax pairs 



L 



n-l 




n > 1, 



(3.8a) 



_ 



'2J / X 



2j 



<Plj 
,V'2j 



U{u, \j 

yW(s,A,) 



, j = 1,2, •■■,7V, 



J = l,2,---,iV; 



l,2,---,iV, 



, j = l,2,^--,iV. 



(3.9a) 

(3.96) 
(3.10a) 

(3.106) 



The spatial part of the nonlinearized Lax pairs and adjoint Lax pairs, i.e. the system (3.9) is a finite 
dimensional system but the time parts of the nonlinearized Lax pairs and adjoint Lax pairs, i.e. the 
systems (3.10) for n > arc all systems of evolution equations in 1 + 1 dimensions. In the following, 
wc shall verify that the system (3.9) is a Liouville integrable (sec Ref. [21]) Hamiltonian system 
and that under the control of (3.9), the systems (3.10) are also Liouville integrable Hamiltonian 
systems. 

The system (3.9) may be cast into the following Hamiltonian form 



dH 



dH 



1,2. 



(3.11a) 



where the Hamiltonian function 



1 2 

H = --{< ^1,^1 >-< $2,*2 >) 

1 2 

- -(< $1,*2 > + < $2,^'l >) - < A$l,*2 > + < A^2,'^l > ■ 



(3.11b) 
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Let us now construct integrals of motion for (3.11). An obvious equality {V^)x = [U,V^] leads to 

Thus F is a generating function of integrals of motion for (3.11). Since F = J2n>o^n^~"', we 
obtain the following expressions 

n 

Fn ^ ^j i,^i^n—i ~t~ ^i^n—i ^i^n— «)• 
i=0 

Further by (3.8), we have 

Fo = - 1, Fi = 2(< $1,^2 > - < $2,*i >), (3.12a) 

n-l 

Fn = {CjCn-i + ajan-i - hK-i) - "^bobn 



i=l 

n-l 



[(< A'-^^i,^i > - < ^^-1$2,^2 >)(< > - < ^"-*-^#2,^2 >) 

i=l 

+ (< *2 > + < A'-^^2, *1 >) (< *2 > + < ^""*"^$2, *1 >) 

- (< A'-^^i, ^2> - < A'-^<^2, ^1>){< A"-^"^^>1, ^-2 > - < ^"-*-^$2, ^1 >) 

+ 2(< ^2 > - < A"-^$2, *i > ), n>2. (3.12b) 

Here F„, n > 0, arc all polynomials of AN dependent variables (pij, tpij. Note that we have (2.10). 
With the same deduction, we find that F = ^trF^ is also a generating function of integrals of 
motion for (3.10). Moreover, we have 

i^.S:f-(-ivAv)MvA:V))^. (3.i3) 



-9$i'5^>2^ V ^ 5$i ^' ^ 5^2 

tr(y-^F) = try r ^^+jAx-^y4r( ^^^^nx-^ 



X-{i+3) 



i>0 
J>1 



2^[q^^-^*i + (Si - 6i)^''"'*2]A-('+^'). (3.14a) 



i>0 



Similarly we can obtain 



tr(F^l^) = 2^[(Si + 6i)^^-'^'i - Q^^-'*2]A-(^+^). (3.14b) 

^ i>0 



The equalities (3.13) and (3.14) lead to 
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The similar deduction can give rise to 



!/(">( |M.„>0. 



At this stage, the systems (3.10) may be readily rewritten as 



d{^Fn+i] 



1 = 1,2. 



(3.15) 



Associated with the symplectic structure on R^^ 



we can generate the Poisson bracket 



{P,Q} = u'^{IdQ,IdP), 



(3.16) 



where IdR denotes a Hamiltonian vector field of a smooth function R on R , defined by IdRA uj = 
—dR, with _1 being the left interior product. Since -F is a generating function of integrals of motion 
for (3.10), we obtain 

(3.17) 



{Fm+i,Fn+i} = 2—Fm+i = 0, m, n > 0. 

This elucidates that {Fn}'^=Q constitutes a Poisson algebra with regard to (3.16). Evidently we 
have 



5F„ 



dF„ 



'E>i='I'2=0 



-2A"-'*i, n > 1. 



(3.18) 



'I'l='I'2=0 



Therefore there must exist one region n C i?4iv ^Yiich the 2N 1-forms dFi,---, dF2N are every 
linearly independent since the Vandermonde determinant V{\i, A2, • • • , \n) is nonzero. In this way 
wc have shown that the spatial part of the nonlinearized Lax pairs and adjoint Lax pairs and the 
time parts of the nonlinearized Lax pairs and adjoint Lax pairs under the control of the spatial 
part are all finite dimensional integrablc Hamiltonian systems in the Liouville sense (sec Rcf. [21]). 



When lim 



\x\ — > + CXD 



^f, = hm 



|x|- 



$j = 0, i = 1, 2 (at this stage, 61 = |Fi = since Fix = 0) 



and $j(a;,t„), ^'j(x,t„) solve simultaneously (3.11) and (3.15), then 



(- < > + < $2,*2 >,- < *1,*2 > - < $2,*1 >) 



is a solution to the nth Dirac system ut^ = Kn- This allows us to conclude that the nth Dirac 
system = has the following involutive representation of solutions 

' q{x,tn) = -< 5&5^„+,$i(0,0),5^5^^^^*i(0,0) > 

' r{x,tn) = - < gfj9'Fl^^<^i{0,0),9^H9Fl+,'^2{0,0) > 

- < 5&5^„+,^2(0,0),5^<7^^^^^i(0,0) >, (3.196) 

with gfj, +1 being the Hamiltonian phase flows (see Ref.[21]) associated with the Hamiltonian 
functions H, F„_|_i, but $^(0,0), 'l'i(0, 0), i = 1,2, being any fixed constant vector of . This 
kind of involutive representation of solutions to the Dirac systems also show, to some extent, the 
characteristic of integrability of the Dirac systems in view of the Liouville integrability of the flows 

4. Conclusions and remarks 

We introduced a kind of symmetry constraint on the Dirac integrable systems. Moreover 
we exhibited an explicit Poisson algebra on the symplectic manifold {R^^ ^uj"^) and further, an 
involutive representation of solutions to the Dirac systems. The proposed symmetry constraint 
includes the eigenfunctions and the adjoint eigenfunctions and is nonlinear with respect to them. 
The corresponding reductions of the Lax pairs and adjoint Lax pairs and the manipulation with 
the integrability of the constrained systems constitute a binary nonlinearization problem. 

We remark that when the zero boundary conditions on the eigenfunctions and the adjoint 
eigenfunctions are not imposed, the spatial part of the nonlincarizcd Lax pairs and adjoint Lax 
pairs is invariant and thus, still a Hamiltonian system with the original Hamiltonian function. 
However, the time parts of the nonlinearized Lax pairs and adjoint Lax pairs will vary with the 
boundary conditions. When the zero boundary conditions are not satisfied, they become very 
complicated Hamiltonian systems. Nevertheless, their Hamiltonian functions are some polynomials 
of F„, n > 1, which may be determined by an obtained equality + — 6^ = 1. Furthermore we 
can consider higher order constraints 

JGn = jJ2'^E/-^, (n>l). 
i=i 

This kind of symmetry constraints is somewhat diflFerent from the Bargmann constraint because 
the Gn involve the diflFerential of u with respect to x. To discuss them, we need to introduce new 
dependent variables, i.e. the so-called Jacobi-Ostrogradsky coordinates. These problems are left to 
a further investigation. 

There exist also certain relations between the nonlinearization technique and the binary non- 
linearization technique. For example, if we take the possible reduction ^'i = —$2, ^2 = ^1, the 
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involutive system (3.12) simplifies to 



Fi =2(< $l,$l > + < $2,^2 >), 

n-1 

Fn=Y^ [(< > - < >)(< ^"-^-^^-l.^l > - < A''-'-^^2,^2 >) 

1=1 

- (< > + < A'-^^2,^2 >)(< > + < ^"-'-^$2,^2 >) 

+4 < ^^-i$i,$2 >< ^"-*-^$i,$2 >] +2(< > + < ^"-^$2,^>2 >), n > 2. 

This is an involutive system on the symplcctic manifold {B?^ ^uP' = d^i A ^$2), and it may be 
generated by the nonlinearization technique. 

We should note that the idea of binary nonlinearization is quite broad; it can be applied to 
other integrable systems (see, for example, Ref . [22] ) . Therefore a large class of finite dimensional 
Liouville integrable Hamiltonian systems may be raised by means of our binary nonlinearization 
technique and the involutive representation of solutions can exhibit the integrability by quadratures 
for the considered systems. 
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